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Abstract 

Using the virial theorem of the heavy-quark effective theory, we show 
that the mixing of the operator for the heavy-quark kinetic energy with the 
identity operator is forbidden at the one-loop order by Lorentz invariance. 
This explains why such a mixing was not observed in several one-loop calcu- 
lations using regularization schemes with a Lorentz-invariant UV regulator, 
and why no UV renormalon singularity was found in the matrix elements 
of the kinetic operator in the bubble approximation (the "invisible renor- 
malon"). On the other hand, we show that the mixing is not protected in 
general by any symmetry, and it indeed occurs at the two-loop order. This 
implies that the parameter of the heavy-quark effective theory is not 
directly a physical quantity, but requires a non-perturbative subtraction. 
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1 Introduction 



The physics of hadrons containing a heavy quark simplifies greatly in the limit 
where the heavy-quark mass mg is taken to infinity. Then new symmetries of 
the strong interactions arise, which relate the long-distance properties of many 
observables to a small number of hadronic matrix elements [jlj]. A systematic 
expansion around the heavy-quark limit is provided by the heavy-quark effective 
theory (HQET), in which a heavy quark inside a hadron moving with four- velocity 
v is described by a velocity-dependent field h v (x) subject to the constraint fh v — 
h v . The field h v (x) is related to the original heavy-quark field Q{x) by a phase 
redefinition, which removes the large "mechanical" part itlqv of the heavy-quark 
momentum arising from the motion of the heavy hadron. Thus, the effective 
field carries the "residual momentum" k = pq — mqv, which characterizes the 
interactions of the heavy quark with gluons. The effective Lagrangian of the 
HQET is!-! 

C cS = h v iv-Dh v + - h v (tD ± ) 2 h v + h v a^GTK + 0(l/m 2 Q ) , (1) 

ZrriQ 4m q 

where = <9 M — ig s A^ is the gauge-covariant derivative, and Dj_ = — 
(v ■ D) v M contains its components orthogonal to the velocity. The origin of the 
operators arising at order l/rriQ is most transparent in the rest frame of the 
heavy hadron: the first operator corresponds to the kinetic energy resulting from 
the motion of the heavy quark inside the hadron (in the rest frame, (iD±) 2 
is the operator for — k 2 ), whereas the second operator describes the magnetic 
interaction of the heavy-quark spin with the gluon field. The Wilson coefficient 
Cmag results from short- distance effects and depends logarithmically on the scale 
at which the chromo-magnetic operator is renormalized. As a consequence of 
the so-called reparametrization invariance of the HQET (an invariance under 
infinitesimal changes of the velocity), the kinetic operator is not multiplicatively 
renormalized |5|, ||. 

The matrix elements of the two dimension-five operators in ([!]) play a most 
significant role in many applications of the HQET. They appear, for instance, in 
the spectroscopy of heavy hadrons and in the description of inclusive weak decays 
[0. For the hadronic matrix elements of the kinetic operator, in particular, one 
defines hadronic parameters Af byQ 0] 

(H(v)\h v (iD ± ) 2 h v \H{v)) = \?, (2) 

where we use a mass-independent normalization of states such that (H(v)\ h v h v \H(v)) = 
1. Spectroscopic relations may be used to extract the difference in the values of 
Af for two hadron states. For the ground-state heavy mesons and baryons, for 

1 Another common notation is to define fi^ H = —X^. 
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instance, we obtain 

~, ™ \ bar \ mos 

^^ = 1 + ^Z^- + (l/ml), (3) 
tub — run ImovriB 

where 0(1 /m? H ) means terms suppressed by three powers of the large hadron 
masses, and Tub = \ (tub + 3m#.) and % = h(mjj + 3rriD*) denote the spin- 
averaged meson masses, defined such that they do not receive a contribution from 
the chromo-magnetic interaction. The value of a single parameter cannot be 
determined from spectroscopy, since it always appears in combination with the 
heavy-quark mass, which by itself is not a physical parameter. This observa- 
tion poses the questions about the status of the heavy-quark kinetic energy as a 
physical parameter. 

In this letter, we study the mixing of the kinetic operator h v (iD±) 2 h v with 
the "identity operator" h v h v under ultraviolet (UV) renormalization. By naive 
dimensional analysis, such a mixing is expected to occur in regularization schemes 
with a dimensionful cutoff parameter A, since the two operators have the same 
quantum numbers but different dimension. This would have important implica- 
tions for phenomenology, as it leads to an additive contribution to the parameter 
Af of the form A 2 C[a s (A)], which must be subtracted in order to define a renor- 
malized ("physical") parameter that is independent of the UV regulator. The 
coefficient C[a s (A)] can be calculated order by order in an expansion in the small 
coupling constant a s (X), and it appears at first sight that the quadratically di- 
vergent term could be subtracted using perturbation theory. This impression is 
erroneous, however, because C may contain non-perturbative contributions of the 
form exp[— 8ir / p a s (\)] = (Aq CD /A) 2 , which cannot be controlled in perturbation 
theory 0. Such terms can contribute an amount of order Aq CD to the parameter 
Af\ which is of the same order as the renormalized parameter itself. Hence, if 
the kinetic operator mixes with the identity operator, it is necessary that the 
quadratically divergent contribution to Af^ be subtracted in a non-perturbative 
way, and hence the heavy-quark kinetic energy by itself is not directly a physical 
quantity. On the other hand, since the quadratic divergence is an UV effect and 
therefore insensitive to the nature of the external states, it follows that the dif- 
ference in the values of Af^ for two hadrons is a physical quantity, and relations 
such as (H) are meaningful. 

The issue of power divergences of matrix elements is intimately related to 
that of UV renormalons |J. The purpose of the heavy-quark expansion is to 
disentangle the short- distance physics characterized by the large mass scale mg 
from the long-distance physics characterized by the typical scale of the momenta 
exchanged between the heavy quark and light degrees of freedom. This is achieved 
by introducing a factorization scale A such that Aq C d "C A <C mq. Contributions 
from momenta above A are controllable in perturbation theory and attributed to 
Wilson coefficients, whereas contributions from momenta below A are contained 
in the matrix elements of the operators in the HQET. If this program is performed 
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with a "hard" factorization scale, these matrix elements diverge, for dimensional 
reasons, as powers of the UV cutoff. For practical reasons, however, one usually 
calculates the Wilson coefficients using dimensional regularization. In this case 
power divergences do not appear, and the factorization scale A is replaced by a 
renormalization scale \i. It is then unavoidable that the Wilson coefficients receive 
contributions from momenta below /i (so-called infrared (IR) renormalons) , and 
the operator matrix elements receive contributions from momenta above fi (so- 
called UV renormalons). These contributions lead to a factorial growth of the 
coefficients in the perturbative expansion of the Wilson coefficients and matrix 
elements. The corresponding perturbation series are divergent and not Borel 
summable; they must be truncated close to the minimal term. Therefore, the 
presence of renormalons leads to intrinsic ambiguities in the definition of the 
Wilson coefficients and matrix elements, which only cancel when all short- and 
long-distance contributions are combined in the heavy-quark expansion [[TTJ-[[nj. 
It is generally believed that from the degree of divergence of the matrix elements 
of an operator in the HQET one can deduce the position of the corresponding 
UV renormalon singularities in the Borel transform of these matrix elements with 
respect to the inverse coupling constant; a power divergence proportional to A™ 
corresponds to a singularity at the position u = n/2 (in a certain normalization 
PI ) on the positive real axis in the Borel plane. For the case of the self-energy of 
a heavy quark, for instance, the correspondence between a linear UV divergence 
and a renormalon pole at u = \ has been established in [T2|, [T3[ . 



The question whether there is a mixing of the kinetic energy with the identity 
operator, and whether there exists a corresponding renormalon singularity, has 
been addressed previously by several authors, with seemingly controversial con- 
clusions. At the one-loop order, such a mixing has indeed been observed when the 
HQET regularized on a space-time lattice f|. Likewise, a "physical" definition of 
a parameter Af^(A) has been suggested, which absorbs certain 0(a s ) corrections 
appearing in the zero-recoil sum rules for heavy-quark transitions fl7| . This defi- 
nition is such that dAf(A)/dA 2 oc a s (A), indicating again a one- loop mixing of the 
kinetic energy with the identity operator. On the other hand, this mixing has not 
been observed at the one-loop order in two Lorentz-invariant cutoff regularization 



schemes [18], which use a Pauli-Villars regulator or a cutoff on the virtuality of 



the gluon in one-loop Feynman diagrams fll9fl . This observation appeared as a 
puzzle, because there seemed to be no obvious reason why the mixing should not 
occur at the one-loop order, and indeed individual diagrams were quadratically 
divergent. Yet, the vanishing of the operator mixing at the one- loop order has 
found its counterpart in the analysis of UV renormalons. As mentioned above, 
a quadratic divergence would correspond to a renormalon singularity at u — 1 
in the Borel plane. Even though we are far from being able to derive an exact 
expression for the Borel transform of a perturbation series, explicit results can 
be obtained in the so-called "bubble approximation", in which all terms corre- 
sponding to an arbitrary number of self-energy insertions on a gluon propagator 
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are kept. Surprisingly, it was found that in this approximation the renormalon 
singularity at u = 1 in the matrix elements of the kinetic operator is absent . 



The absence of the related IR renormalon in the pole mass of a heavy quark was 
noted previously in ||12|| . This "missing" singularity has been called the "invisible 
renormalon" . Three possible explanations for this puzzle have been suggested: 

1. The vanishing of the operator mixing at the one-loop order is accidental 
and happens only in some peculiar regularization schemes. However, this 
does not explain why the renormalon singularity at u — 1 is "invisible" in 
the bubble approximation. 

2. The vanishing of the mixing is an artefact of one-loop perturbation theory. 
The mixing appears in higher orders, and the corresponding renormalon 
singularity appears when one goes beyond the bubble approximation. 

3. There is a symmetry which prevents the mixing between the kinetic en- 
ergy and the identity operator, but this symmetry is broken by the lattice 
regularization. Candidates of such symmetries are Lorentz (or rotational) 
invariance and the reparametrization invariance of the HQET || ||. In 
fact, it has been conjectured that Lorentz invariance may protect the ma- 
trix elements of the kinetic energy against an UV renormalon singularity at 
u = 1 P4]| . However, later it was argued that there is no obvious connec- 
tion between Lorentz or reparametrization invariance and the absence of a 
quadratic divergence (or the corresponding UV renormalon) in the matrix 
elements of the kinetic energy [0 . 

In this paper we shall resolve this puzzle. We will prove that the mixing 
between the kinetic energy and the identity operator is forbidden at the one- 
loop order in all regularization schemes with a Lorentz-invariant UV cutoff. For 
our arguments to hold, it is essential that the regulator be introduced in a fully 
Lorentz-invariant way; in particular, it is not sufficient to choose a regularization 
scheme that respects only rotational invariance in the rest frame of the heavy 
hadron. The fact that a mixing at the one-loop order was observed in || and 
H] is a consequence of the explicit breaking of Lorentz invariance by the regu- 



larization schemes adopted in these calculations. On the other hand, in general 
there is no symmetry that protects the matrix elements of the kinetic operator 
from quadratic divergences, and we show explicitly that a mixing with the iden- 
tity operator occurs from the two-loop order on. In the renormalon language, 
this suggest that there is a singularity at u = 1 in the Borel plane, which would 
appear when one goes beyond the bubble approximation. Hence, each of the 
above-mentioned suggestions is only partly correct, and the puzzle of the "invis- 
ible renormalon" finds its explanation in an "accidental" space-time symmetry, 
which is only operative at the one-loop order. 
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Figure 1: One- loop diagrams that could contribute to the mixing of the kinetic 
energy with the identity operator. The kinetic operator is represented by a grey 
square; heavy-quark propagators are drawn as double lines. 



2 Operator Mixing 

As mentioned above, in regularization schemes with a dimensionful cutoff A, the 
kinetic operator can mix with the identity operator under UV renormalization. 
We define a renormalization constant Z D2 _+\ by 



h„ (iD i ) h, 



^ = Z D2 ^i h v h v + ... , (4) 



where the subscript "bare" indicates the bare operator, whose matrix elements 
depend on the UV regulator, and the ellipses represent contributions from op- 
erators of higher dimension. Such a mixing leads to an additive contribution to 
the matrix elements of the kinetic operator which, for dimensional reasons, is 
quadratically divergent in the cutoff A. It follows that 

1 Af(A) = A Zfla _ +1 . ( 5 ) 



dA 2 1 v ' dX 

Our goal is to perform a perturbative calculation of Z D2 _^i in a regularization 
scheme where the cutoff preserves Lorentz invariance. By evaluating the UV di- 
vergences of the diagrams shown in Fig. [I] (supplemented by the wave-function 
renormalization of the quark fields), it has been shown in |TB[ that Z D 2^\ vanishes 
at the one-loop order in two such regularization schemes. However, individual 
diagrams have a quadratic divergence, and the reason for the cancellations oc- 
curring in the sum of all diagrams remained unclear. 

To extend the direct calculation of Z D2 _>i to the next order would involve the 
evaluation of a large number of two-loop diagrams. In this letter, we suggest a 
more efficient way to perform this calculation, which in addition will provide an 
explanation for the vanishing of Z D %_^ at the one-loop order. To this end, we use 
the virial theorem of the HQET, which establishes a relation between the kinetic 
energy and a matrix element of an operator containing the gluon field-strength 
tensor between states moving at different velocities f7|, |2(J . For the ground-state 
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Figure 2: One-loop diagrams that could contribute to the mixing of the chromo- 
electric operator with lower- dimensional operators. The velocity- changing 
chromo-electric operator is represented by a white square. 



pseudoscalar mesonsQ and Aq baryons, the virial theorem states that 

(H(v')\ h v ,ig 8 G^h v \H(v)) = OV* - vV) |-^- + 0(v ■ v' - 1) j . (6) 
This is a non-perturbative relation between hadronic matrix elements, which is 



preserved by renormalization [21]. In the rest frame of one of the hadrons, only 
the electric field has a non- vanishing matrix element, and we shall therefore refer 
to the operator on the left-hand side as the chromo-electric operator. Note that 
(H) is not an operator identity, but rather a relation between ground-state matrix 
elements; the corresponding operator identity would contain operators of higher 
spin, which have vanishing ground-state matrix elements. However, ([]) is an 
operator identity as far as the quadratic UV divergences are concerned. 

Similar to the kinetic operator, also the chromo-electric operator can mix with 
a dimension-three operator. Since the Feynman rules of the HQET do not involve 
Dirac matrices, there is only one candidate for such an operator, and we define 

h v ,ig s G^h v = Z E ^i(v ■ v') (v»v' u - v u v'») h v ,h v + . . . . (7) 

Although there exist other dimension-three operators with the right quantum 
numbers, such as h v 'ia^ u h v and h v '[(v — v'Y^ u — {y — v') u ^] h v , they cannot be 
induced by quantum effects. The virial theorem (^]) implies a relation between the 
renormalization factors Z E ^i and Z D2 _^i in the limit of equal velocities (v -v' = 1), 
which reads 

Z D2 ^ = -3Z E ^(1). (8) 

Thus, if we understand the mixing of the chromo-electric operator in (0), we also 
understand the mixing of the kinetic energy with the identity operator. As we 
shall see, calculating Z D2 ^i by means of the relation (||) is much more efficient 
than a direct calculation. 

Since the gluon field-strength is an antisymmetric Lorentz tensor, the only 
possible dimension-three operator appearing on the right-hand side of (^) contains 

2 The same relation holds for the ground-state vector mesons, when an averaging over the 
two transverse polarization states is implied. 
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both velocities, v and v'. Consider now the diagrams shown in Fig. 0, which could 
potentially give rise to operator mixing at the one-loop order. Let us assume that 
we regulate the UV divergences of these diagrams in a Lorentz-invariant way, for 
instance, by introducing a Pauli-Villars cutoff. Since one leg of the gluon must be 
attached to the chromo-electric operator, there is no way to obtain a dependence 
on both heavy-quark velocities; the first diagram can only depend on v, while 
the second one can only depend on v' . Hence, because of Lorentz invariance 
the renormalization constant Z E _^i must vanish at the one-loop order, and so 
does Z D 2^\. The relation (f|) implied by the virial theorem thus provides for a 
very simple explanation of the apparently accidental cancellations between the 
quadratic divergences of the one-loop diagrams in Fig. [I], observed in |18| . 



On the other hand, the above argument is clearly linked to the one-loop order. 
Let us, therefore, study the problem at the level of two loops. Then there are 
four diagrams that could potentially lead to a dependence on both heavy-quark 
velocities. They are shown in Fig. [3]. (Notice that this is a vast simplification over 
the direct calculation of the two- loop mixing of the kinetic operator.) We need to 
calculate the UV divergences of these diagrams in a regularization scheme with a 
dimensionful cutoff. However, we have to be careful not to spoil Lorentz or gauge 
invariance when introducing the cutoff. For the case at hand, this can be done 
using dispersion relations. Let us consider the matrix element of the chromo- 
electric operator evaluated between quark states with residual momenta k and 
k! . (The nature of the external states is irrelevant, since we are only interested 
in the UV behaviour of the diagrams.) This matrix element has discontinuities 
if either v ■ k or v' ■ k' are positive, since then physical intermediate states can 
be excited. The Feynman amplitudes corresponding to the diagrams in Fig. [3] 
satisfy a dispersion representation of the form[] 

A=[dtO /do/, : r, (9) 

J J (co-v- k-ie){u'-v'- k'-ie) ' K1 

where p(u, to') is the double spectral density. The variables to and to' have a 
concrete physical meaning as the energies of the intermediate states, measured 
with respect to the heavy-quark mass. A mixing with the identity operator 
would show up as a quadratic divergence of the dispersion integral arising from 
the region of large energies. Introducing an UV cutoff in the dispersion integral 
regulates these divergences without spoiling any of the symmetries of the theory; 
in particular, Lorentz and gauge invariance are preserved. 

Below, we shall investigate two different prescriptions how to regulate the 
dispersion integral. The most direct way is to introduce a hard UV cutoff, i.e. 

a^= LL^. (io) 



3 Once the double integral is regulated in the UV region, there are no subtraction terms in 
the dispersion relation. 
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v,k v , k v v v v v v 

Figure 3: Two- loop diagrams that could contribute to the mixing of the chromo- 
electric operator with lower-dimensional operators. Not shown are two copies of 
the first two diagrams with the "outer" gluon attached to the other heavy-quark 
line. 

As we are interested in studying the UV behaviour of the operator matrix element, 
we have set the external momenta to zero (there are no IR divergences). For 
dimensional reasons, the result will then be quadratic in the cutoff A. Another 
scheme is to perform a double Borel transformation of the original integral in the 
external variables v ■ k and v' • k'. This leads to 

A Borel = _J_ f duJ f dc y p(ll , J) e -»/M-u>/M> (U) 

MM' J J v ' 

corresponding to a "soft" exponential cutoff. It is most natural to set the Borel 
parameters equal: M = M' = A, in which case the result will be proportional to 
A 2 . 

To evaluate the two-loop diagrams in Fig. [3], we work in the Feynman gauge; 
the sum of all diagrams is gauge independent. Using the Feynman rules of the 
HQET, it is then straightforward to see that the first two diagrams vanish. The 
third diagram is easily calculated, since it factorizes into two one-loop integrals. 
Using dimensional regularization, with d being the number of space-time dimen- 
sions, we obtain 

-2C F C A -^-7 r2 (3 - d) Y 2 (d/2 - 1) (-2v ■ k) d ~ 3 (-2v' ■ k') d ~ z 

(471")" 

x(v»v' v -v v v'»)u v ,(k')u v (k). (12) 

Here Cf = \{N 2 — 1)/N C and Ca = N c are the eigenvalues of the quadratic 
Casimir operator in the fundamental and the adjoint representations, and u v (k) 
is a heavy-quark spinor. The external momenta k and k! act as IR regulators. 
The result is real if v ■ k < and v ' ■ k' < 0. However, as expected there are 
discontinuities for positive values of these variables. The corresponding double 
spectral density is finite in four dimensions and given by 

a 2 

~ C F C A tt^ 9(u) 9(u') ujuj' (vN v - v v v'») u v ,{k') uJk) . (13) 
2ti 

The calculation of the last diagram is more involved, since this is a genuine two- 
loop diagram with a non-trivial dependence on the variable v ■ v'. We shall first 
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discuss the result obtained in the scheme based on using a double Borel transfor- 
mation. For this case, efficient techniques to calculate two-loop integrals involving 
two heavy quarks moving at different velocities have been developed in [22]. In 
the case of equal Borel parameters (M = M' = A) we find, after combining the 
results for the two diagrams, that the renormalization factor defined in (0) is 
given by 



o: 2 



Z%$(y ■ v') = -C F C A ^ F(v ■ v') A 2 , (14) 



where 



I' Or) = \L 2 (-w + )-L 2 (-w^)]+^p^ hi[2(«- l)J - — — . (15) 



{w 2 - If/ 2 



itr — 1 w 



Here w± — w ± y/w 2 — 1, and L 2 (x) = — J a ^ln(l — y) is the dilogarithm. 
Evaluating this expression in the equal- velocity limit, using that F{1) = |(1 + 
2 In 2), we obtain from (|5p the two- loop renormalization constant for the mixing 
of the kinetic energy with the identity operator. The result is 



Z^tl = CfC a || (1 + 2 In 2) A 2 . (16) 

The precise numerical coefficient of the quadratic divergence depends on the 
way in which the UV cutoff A is introduced; below we will derive the result for 
another scheme. However, the fact that we have found a regularization scheme 
that preserves Lorentz and gauge invar iance, and in which the kinetic energy 
mixes with the identity operator at the two-loop order, proves that such a mixing 
is present in general. Beyond the one-loop order, it is not forbidden by any 
symmetry. 

Let us now elaborate on the result (0). We have calculated the generalization 
of this relation to the case of unequal Borel parameters. Then 
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Z^i = C F C A ^ MM' G(M/M') , (17) 
where 

, 2(2 - 3x + 2x 2 ) , M . 2a; 2 (l-3a;-2x 2 ) , 1 - Qx + x 2 /io . 
G(x) = — Mn l + x + K — ^ lna; ^— . 18 

w x v ' (1 + x) 3 (l + x) 2 y ' 

The function G(x) is symmetric under the exchange x «-> 1/x and obeys G(l) = 
1 + 2 In 2. Expression (0) is useful, since it allows us to derive the double spec- 
tral density in the dispersion integral @. It is obtained by performing another 
double Borel transformation, now in the variables —1/M and —1/M', with Borel 



parameters 1/lu and 1/lu' ||23||. The calculation is straightforward and leads to 
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where 

p(u, (J) = 9(oj) 9{u') Ljoj' + 5c; 2 , + 2(3uu' - u 2 - u' 2 ) In ^1 - ^ | , (20) 
with u < = min(u,u') and u > = max(w,w'). With the explicit result for the 



double spectral density at hand, we can now evaluate, according to fllOf) , the 



renormalization factor in the scheme with a hard UV cutoff. The result is 



a 2 /7r 2 



Z oLi = C F C A ^- 2 {--l)\ 2 . (21) 
As expected, the numerical coefficient in front of the quadratic divergence is 



different from the case of the Borel regulator in (16), but the general structure 
of the result is the same in both regularization schemes. 

We emphasize that in deriving that there is no mixing at the one-loop order 
it is crucial that the UV cutoff be introduced in a fully Lorentz-invariant way; in 
particular, for our arguments to hold it is not sufficient to choose a regularization 
scheme that respects only rotational invariance in the rest frame of the heavy 
hadron. This is relevant when comparing our results to some previous calcula- 
tions, in which a non- vanishing contribution to Z D - 2 ^i was obtained at 0(a s ). In 
the case of the lattice regularization, the one-loop mixing of the kinetic energy 
with the identity operator arises because the lattice breaks Lorentz invariance 
(i.e., Euclidean rotational invariance). The result for the renormalization factor 
Z D 2^\ in this scheme is || 



^n^Li — ~~ fciatt Cf — A 2 , (22) 

71 



where A = n/a is the UV cutoff (a is the lattice spacing), and the numerical 
coefficient k\ att is given by 



_L_ 

^latt n t 
ATX 6 



1 3 

d 3 k ; , A = Y (1 - cos hi) . (23) 

''1 + A) 2 -1 



i=l 



The appearance of the quantity A shows that the lattice regularization breaks 
both Lorentz and rotational invariance. Things are slightly more intricate with 
the scheme adopted in |17]], however. The definition for the parameter A^(A) 
suggested there is equivalent to introducing a cutoff on the spatial momentum 
(|k| < A) in the rest frame of the hadron containing the heavy quarkfj Although 
working in the rest frame is quite natural when dealing with heavy quarks, this 
prescription breaks Lorentz invariance in that it specifies a particular reference 
frame, thereby introducing a dependence on the velocity v through the regular- 
ization procedure. Then the arguments that led us to conclude that the one-loop 

4 The author is grateful to M. Shifman for pointing out this fact, and for suggesting the 



derivation of (26) along the lines discussed here. 
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diagrams in Fig. ^ have no quadratic divergences are no longer valid. The second 
diagram, which explicitly depends on v', now acquires an implicit dependence on 
v through the regularization, and thus it can contribute to Z E _^i in (|7]). After a 
straightforward calculation, we find that 

rv 

Z E ^(vv') = C F7 ^H(vv')X\ (24) 

where 

H(w) = w ' r ^\ r (w)= - I . ln(w + Vw^l) . (25) 
w z — 1 y/w 2 — 1 v ' 

The same result is obtained if the cutoff is introduced in the rest frame of the 
final-state hadron. Using that H(l) = |, we obtain from @ 

Z D2 ^ = -C F ^ A 2 , (26) 

7T 

which agrees with the result derived in |I7|]. The same result can be recovered 
from (]22D by taking the "continuum limit" A — > ~ k 2 and imposing the cutoff on 
|k| (i.e. by integrating over a sphere rather than a cube), in which case rotational 
invariance is restored, and fci at t — > 1- Since the lattice version of the HQET is 
formulated for heavy quarks at rest, Lorentz invariance is broken even in the 
"continuum limit" . 



3 Conclusions 

Using the virial theorem of the heavy-quark effective theory, we have shown that 
the mixing of the operator for the heavy-quark kinetic energy with the identity 
operator is forbidden at the one-loop order by Lorentz invariance. This resolves 



the puzzle of the "invisible renormalon" raised in [jig] , which consisted in the 
observation of an apparently accidental vanishing of the mixing at the one-loop 
order in several Lorentz-invariant cutoff regularization schemes, and the asso- 
ciated absence of an UV renormalon singularity in the bubble approximation. 
However, the mixing is not protected in general by any symmetry, and we have 
shown explicitly that it does occur at the two-loop order. This means that the pa- 
rameter Af^ of the heavy-quark effective theory is not directly a physical quantity, 
but requires a non-perturbative subtraction. 

Our results for the two- loop renormalization of the kinetic energy imply that, 
in any regularization scheme with a dimensionful cutoff parameter A that pre- 
serves Lorentz and gauge invariance, there is an additive, quadratically divergent 
contribution to the parameter A^(A) of the HQET, such that 

Af(A) = = * C r C A^zr + 0(a 3 s ). (27) 



dA 2 'I V V u^l 47f2 
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For the two regularization schemes considered here, the numerical coefficient 
k takes the values (1 + 2 In 2) « 2.39 and (±tt 2 - 1) « 2.29. Note that this 



renormalizat ion-group equation is different from the one derived in [17], which 



contains a one-loop contribution of the form —Cf a s (X)/rr on the right-hand 
side, in accordance with (|26|). The reason is that the "physical" definition of 
the parameter /J% H (\) = — Af*(A) suggested by these authors is equivalent to 
introducing a cutoff on the spatial momentum in the rest frame of the heavy 
hadron and thus breaks Lorentz invariance. (Although it is perfectly consistent, 
this definition may therefore not be the most convenient one.) The values of the 
kinetic energy defined in different regularization schemes can be related to each 
other using perturbation theory. From our results, it follows that at the one-loop 
order there is a universal relation between the parameter ^ H {X) defined in IT7J 
and Af*(A) defined in any Lorentz-invariant regularization scheme, which reads 

A*.* (A) = "Af (A) + C F ^1 A 2 + 0{a]) . (28) 

7T 

This implies that the lower bound for the kinetic energy of a heavy quark inside 
a meson, fJ^ me8 (X) > 0.36 GeV 2 , which by construction holds for the defini- 
tion proposed in fL7| , becomes weaker when the kinetic energy is regulated in a 
Lorentz-invariant way. Even if the cutoff A is chosen as low as 1 GeV, the extra 
term in ( p8"D amounts to about 0.15 GeV 2 , implying that the lower bound for 
— A™ es (A) is at most 0.2 GeV 2 . This observation is in accordance with arguments 
presented in [EMI, where a similar conclusion has been reached. 



As mentioned in the introduction, the presence of a quadratic divergence en- 
tails that the parameter A^ of the HQET requires a non-perturbative subtraction 
to be defined and as such is not directly a physical parameter. Since in Lorentz- 
invariant regularization schemes the operator mixing occurs only at the next- 
to-leading order, one may hope that in practice the ambiguities related to the 
choice of the subtraction scheme will be small. Nevertheless, one has to be careful 
when comparing the values of Af* obtained using different theoretical methods. 
For instance, QCD sum-rule determinations of Af* for the ground-state heavy 
p3|— 1|27|| and baryons [ES| are performed in a scheme based on dimen- 



mesons 



sional regularization, i.e. without a dimensionful UV cutoff. Strictly speaking, 
the results are then affected by a renormalon ambiguity problem, i.e. there is an 
intrinsic arbitrariness of order Aq CD , which would show up when higher-order 
perturbative corrections to the sum rules would be calculated beyond the bubble 
approximation. Still, these results can consistently be used in connection with 
Wilson coefficients calculated to the same accuracy. Therefore, the sum-rule cal- 
culations may be compared with the results of phenomenological determinations 
of A° ies from a combined analysis of inclusive decay rates and moments of the de- 
cay spectra in decays of beauty and charm mesons p9|-[5I|. On the other hand, 



when Af* is determined using lattice gauge theory, the presence of a dimensionful 
regulator (the lattice spacing) requires a non-perturbative subtraction, and the 
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results depend on how this subtraction is performed. A "physical prescription" 
for such a subtraction has been suggested in . 
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